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1 INTRODUCTION 



ABSTRACT 

In an ongoing effort to study the formation and evolution of galactic bulges we have 
investigated the position of 19 bulges of type SO-Sbc on the fundamental plane (FP). 
We find that bulges, both in B and K band, lie close to but slightly below the FP 
defined by ellipticals and SOs, i.e. are slightly brighter. There are hints that bulges 
of later morphological type are situated below the other bulges in our sample. The 
FP results are consistent with the picture, obtained from our recent analysis of HST- 
colours, that bulges are old, except for the Sbc galaxies. The fact that bulges lie so 
close to the FP of ellipticals and SOs implies that their formation epoch must have 
been similar, at most about 2.5 Gyr earlier, than cluster Es and SOs, and and that the 
surrounding disk does not significantly affect their structure. 

Key words: galaxies: bulges - galaxies: fundamental parameters - galaxies: indi- 
vidual(NGC5326, NGC5389, NGC5422, NGC5443, NGC5475, NGC5587, NGC5689, 
NGC5707, NGC5719, NGC5746, NGC5838, NGC5854, NGC5879, NC5965, NGC6010, 
NGC6504, NGC7332, NGC7457, NGC7537) 
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(1987) and Dressier et al. (1987), the Fundamental Plane 



initial mass functio n (IMF) or even the dark matte r distri- 



(FP) establishes llial Lire effeelive ladius (r e ), the mean ef- 



bution in galaxies (Suzman, Lucey fe Bower 1993 



fective surfac e brightness (E e ) and the central velocity dis- 

persion (oo) are related by 

r e = C ■ cr ■ L e 

The FP is one of the main scaling relations for early- 
type galaxies (e.g. Pahre, Djorgovski & de Carvalho 1998, 
hereafter P98). The fact that, compared to other relations, 
its scatter is so small, has made it a useful distance indicator 
(J0rgensen, Franx & Kjaergaard 1996, hereafter J96). Histor - 
ically, the FP is a refinement of the Faber- Jackson (1976), 
relation between galaxy luminosity and ctq, which includes 
the depende ncy o f luminosity and <ro on surface brightness 
(Kormendy (tl985h). 



Lanzoni & Renzini 1996) 



Ciotti 



Beyond its value as a distance indicator, the FP con- 
tains useful clues on the physics of galactic systems. The 
position of galactic bulges on the FP of elliptical galaxies 
provides information about the differences in stellar popu- 
lations, internal structure and kinematics of both types of 
galaxies. In a series of papers, Prugniel & Simien (1994, 
1996, 1997) have gone into the details of using the FP for 
this purpose. Assuming that the FP is a representation of 



the Virial Theorem a 



— , and that the mass-to- light 



Much work has been devoted to determining the expo- 
nents A and B. The virial theorem predicts A = 2 and 
B — — 1. However, observations show lower values of A 
and higher values of B (Busarello et al. 1997). The value 
B has been well determined to be ~ -0.8 (P98, J96). A, 
usually referred to as the 'tilt' of the FP (P98), has val- 
ues that vary significantly depending on the band employed. 
Several mechanisms are likely to produce the difference be- 



ratio scales as M- oc L 13 , these authors obtain the relation 
between observables <tq k —L 1+i3 . In their picture, this 
relation is tight for old elliptical galaxies, since for them cro 
is a good mass-indicator, and is a smooth function of L, 
probably metallicity-driven. Deviations to the relation can 
be due to the presence of young stellar populations (causing 
a change in see also Schweizer & Seitzer 1992), rota- 
tional support (since in this case oo gives an underestimate 
of the gravitational energy), or differences in structure, for 
example indicated by the shape of the surface brightness 
profile. 



© 0000 RAS 



2 J. Falcon- Barroso, R.F. Peletier, M. Balcells 



One of the first papers addressing the position of bulges 
on the FP is the study of Bender, Burstein & Faber (1992, 
hereafter B92). They present B-band FP properties for a 
large sample of dynamically hot galaxies, which they di- 
vide into several groups, referred to as Es, intermediate 
Es, compact Es, bulges, bright dwarf ellipticals and dwarf 
spheroidals. B92 conclude that bulges lie preferentially be- 
low the FP of ellipticals, suggesting that, on average, bulges 
are younger than ellipticals. Most (17 out of 19) of the bulges 
in B92 are in lenticular galaxies. In a related study using 
a sample of gal axies from the Coma cluster Saglia, Ben- 
der & Dressier (1993) find a small shift in the same direc- 
tion in the location of bulges when compared with Es and 
SOs (0.1 in log(r e ) in B-band). In contrast with these find- 
ings, J96, in their comparison between the FP of 10 clus- 
ters, find a negligible difference between E and SO galaxies 
(A(E - SO) = 0.01 ± 015 in B-band), a result shared by 
Prugniel & Simien ( |l99e| ). 



The FP of bulges is more complicated to study than 
that of elliptical galaxies, since a determination of struc- 
tural parameters requires a decomposition of the light into 
bulges and disks, as well as to take into account dust red- 
dening and extinction. These problems only worsen as we 
move to later types. Perhaps as a result of this, very lit- 
tle information is available about the FP of bulges of types 
later than SO. This paper intends to fill this gap by present- 
ing a study of the FP of bulges for types ranging from SO 
to Sbc. It belongs to a series of papers d iscussing the prop- 



erti es of galactic bulges. The first paper (Balcells & Peletier 
1994Q 

iscusses optical colours and colour gradients of an 
optically selected, complete, sample of ~ 50 bulges of type 
SO-Sbc. The galaxies were required to have an inclination i 
>50° to minimize the amount of disk light which projects 
onto the bulge and ensure that extinction from dust in the 
disk does not dominate the colours measured on one side of 
the bulge. Of this sample, 30 galaxies for which the colour 
profiles were smooth on one side of the bulge (i.e. where the 
colours were not significantly affected by dust extinction) 
were observed in the near-infrared iC-band at higher spatial 
resolution. Andredakis, Peletier & Balcells (1995, hereafter 
APB95) present a photometric decomposition into bulge and 
disk, using Kent's (1986) model-independent method, after 
which Sersic r 1 ^" laws are fitted to the bulge. They find 
that the Sersic law generally fits very well, and that there 
is a good correlation between n and morphological type, 
such that bulges of early-type spiral galaxies follow the r 1//4 
law, while profiles of late type spirals are more nearly expo- 
nential. Optical and near-infra red co lours of bulges are pre- 
sented in Peletier & Balcells (1996). Photometric data for 



this sample, including colour profiles along major and minor 
axis and fits of ellipticity and position angle as a function of 
radius, are presented in Peletier & Balcells (1997, hereafter 
PB97). Central color profiles, based on HST/WFPC2 and 
NICMOS imaging of a su bsamp le of 20 of our galaxies, are 
analysed in Peletier et al. (1999). The latter paper concludes 



that all bulges, except those of Sbc's, are old with an age 
spread of less than 2 Gyr. The present paper makes use of 
our own spectroscopy of the sample to address the position 
of bulges on the FP of elliptical galaxies. A full report of the 
kinematic properties derived from the spectroscopy will be 
given in a forthcoming paper. 



When analyzing bulges, several authors have made the 
assumption that bulges' surface brightness p rofiles are ei- 
ther r 1 / 4 or exponential ( Courteau et al. 1996| ), a dichotomy 
of shapes that would suggest a dichotomy of formation pro- 
cesses. Carollo ( |1999| ) also divides bulges into r 1//4 and expo- 
nential classes in her study of galaxy nuclei. In our work we 
define and characterize a bulge on the basis of its signature 
in the galaxy's ellipticity profile, rather than on the shape 
of the surface brightness profile. We thus avoid having to 
assign each bulge to a given profile class, and indeed we find 
a continuum of profile shapes, which we model with Sersic's 
law with free index n, rather than two separate classes. Be- 
cause our sample includes objects ranging from exponential 
to r 1 ' 4 , our FP study should uncover any structural and 
population signatures that depend on the shape of the sur- 
face brightness profile. 

Section §2 summarises the observations and data reduc- 
tion as well as procedures to derive the different structural 
and kinematic parameters. In §3 we show the FP in the dif- 
ferent bands and compare the results with other authors. 
The Mg2 - <7 relation is analysed in §4. Possible correlations 
between observables and the residuals of the FP are studied 
in §5. Throughout the paper we use a Hubble Constant of 
#0=50 kms- 1 Mpc _1 . 



2 OBSERVATIONS & DATA REDUCTION 
2.1 Observations 

We obtained long-slit spectra of our sample galaxies at the 
4.2m William Herschel Telescope of the Observatorio del 
Roque de los Muchachos at La Palma between 11 and 13 
July 1997. Nineteen galaxies were observed (see Table 1). All 
galaxies were observed between 3660 - 5560 A and between 
8360 - 9170 A using the ISIS blue and red arms, respectively. 
The red arm was equipped with the Tektronix ( 1024 x 1024) 
TEK2 CCD (0.36" per 24 fim pixel) and the R600R grat- 
ing, giving a spectral resolution (FWHM of arc lines) of 
1.74 A (60 km/s). The blue arm was equipped with the Lo- 
ral (2048x2048) LOR2 CCD, (15 /im pixels, 0.22 "/pixel) 
and the R300B grating, giving a spectral resolution of 4.1 
A. We used the 6100 A dichroic for simultaneous red and 
blue arm exposures. Typical exposure times per galaxy were 
1500 sec in both wavelength ranges. The slit width was 1.2", 
matching the seeing at the time of the observations. Arc line 
exposures were taken before and after each target exposure. 
Tungsten continuum lamp exposures were taken with the 
red arm after each target exposure, for fringe calibration. 
Twilight sky exposures were taken every night for flat field- 
ing. We took spectra of spectr al templates of type s A2 to 
M5, from the Lick list of stars ([Worthey et al. 1994|). 



2.2 Spectroscopic Parameters 

2.2.1 Spectral Data Reduction 

The data were reduced following standard steps using the 
IRAF package. Each individual frame was divided by a flat- 
field, taken at the same position, after a bias frame had 
been subtracted. We subtracted the sky using the outer re- 
gions of the 4' long slit, an d removed the cosmic rays us- 
ing the REDm C E package ( |Cardiel 1999|) . The data were 
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Table 1. The Data Set. 
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NOTES: Photometric errors in columns (5) and (7) are 0.05 and 0.08 mag respectively (PB97). For definition of log(cro) an d Mg2,o see 

text. 



kindly made available to us by the authors. A Gaussian fit 
to the line-of-sight velocity distribution (LOSVD) was per- 
formed to obtain mean radial velocities and velocity disper- 
sions. Star templates of spectral types G8III, KHIIb, K3III 
and K5III, gave dispersions with rms differences around 5 
km/s. We took the median value of the 4 determinations 
as the best estimate of the central velocity dispersions of 
our galaxies. Velocity dispersions results are given in col- 
umn (8) of Table |l| Aperture corrections have been applied 
to the central velocity dispersions. We have converted a rect- 
angular aperture (1.1 x 1.2 arcsec) to an 'equivalent circular 
aperture'. Then we have corrected those values to a stan- 
dard aperture, in this case denned to be 1.7 arcsec at the 
distance of Coma (see j0rgensen et al. 1995 for more details) 
The velocity dispersion errors in Table ^ are uncertainty 
estimates derived from simulations using FOURFIT. Sim- 
ulated galaxy spectra were made by broadening a K3III 



calibrated in wavelength, with an RMS error in the calibra- 
tion solution of 0.1 A in the blue and 0.02 A in the red 
arm. Finally the d ata were c alibrated in flux u sing standard 
stars BD+284211 ( pke 1990[ ) and BD+174708 ( |Oke fc Gunn 



19830 



2.2.2 Derivation of kinematic parameters 

Kinematic parameters were derived from data in the region 
of the IR Ca II triplet (8498 A, 8542 A, 8662 A) in the red 
arm spectra. We worked with coadded spectra covering the 
central 3 pixels of each long-slit galaxy spectrum, roughly 1.1 
arcsec, which is a good match to the seeing conditions dur- 
ing the observations. Signal-to-noise ratios per pixel for these 
spectra are listed in Table |l[ We made use of the pr ogram 



FOURFIT developed by van der Marel & Franx ( |1993| ) 
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star with different LOSVD, using Gaussian as well as non- 
Gaussian profiles, and adding white noise to yield S/N of 
25, 33, ,50 and 100 These simulated spectra were analyser! 



I(r) = I e ■ 10 



(-^■([r/r.] 1 /"-!)) 



(1) 



with t hlree star templates of types G8TTT, G9VT and K3TTT. 
For each S/N we define the dispersion error as the rms of all 
the velocity dispersions coming from all templates as well as 
all LOSVDs. The error estimate thus includes contributions 
from S/N, template mismatch and non-Gaussian LOSVD. A 
complete report of the simulations that we performed will be 
given in a forthcoming paper that presents the full LOSVDs 
of the galaxies. 



2.2.3 Derivation of Mg 2 ,o 

The blue arm data was used exclusively to measure the cen- 
tral Mg2 index, here indicated with Mg2,o- We have mea- 
sure d the aperture corrected Mg2,o index (see j0rgensen et 
al. (1995 ). for more details) and converted it to the Lick 
system ( Worthey et al. 1994 ) following the procedure in 
Vazdekis (|l997|)! 



2.2.4 Derivation of distances 

Distances to the galaxies wer e derived from optical reces - 
sion velocities given in RC3 ( ]de Vaucouleurs et al. 199l|). 
corrected to the centroid of the Local Group ( [Karachent- 
3CV fc Makarov 1996| ), and assuming a uniform Hubble flow. 
Distance errors are therefore significant, due to the Virgo- 
centric flow and the possible influence of the Great Atractor. 
We have taken 150 km/s as the error in the radial velocity. 
We ha ve taken this value from the work of Aaronson et al. 
( |1982| ) who studied the influence of bulk motions on a sam- 



ple of spiral galaxies at distances comparable to our sample. 



2.3 Photometric Parameters 

2. 3. 1 Derivation of structural parameters in K-Band 

The derivation of the structural parameters in K-band was 
done as follows. The K-haxid effective radii in arcsec were 
obtained from APB95. These r e therefore were derived from 
Sersic fits to if-band bulge profiles after an ellipticity- 
based bulge-disk decomposition (see APB95 for details). The 
APB95 major-axis r e values were scaled to the geometric 
mean radius r e = ^ 

Te, minor ' fe, major for direct comparison 
with other authors (B92, J96, P98) using ellipticity profiles 
from the K-band ellipse fits published in PB97. 



2.3.2 Derivation of structural parameters in B-Band 

For the B-band FP, our goal is to minimize the effects of dust 
extinction on the photometry. An r Ei s derived from a direct 
ellipse fit to the B-band image is likely to be affected by the 
copious dust pres ent on one side of t he galaxy, and also in 
the nuclear area (Peletier et al. 1999). We opt therefore to 
derive the dust-corrected _B-band effective radius r e ,s from 
the effective radius in K and the colour gradient determined 
on the minor axis. The conversion was done in the following 
way: 

We assume that our profile is of the form: 



with b n = 0.868 ■ n - 0.142 (Caon, Cappacioli & D'Onofrio 



1993j ). Differentiating B - K = -2.5 ■ log[I B (r)/I K {r)] + est 
with respect to log(r) gives: 

d{B- K) _ r d(B — K) 



d log(r) log(e) 



dr 



(2) 



Working this expression out, we get the relation be- 
tween the effective radii in the two bands as: 



r e ,B 



X7(B-K) 
S(n) 



' r e ,K 



(3) 



where V(B — K) = d(B — _K")/dlog(r) is the colour 
gradient taken from PB97 and S(n) = (2.5 ■ b n )/(n ■ log(e)). 
The ratio r e ,s/r e ,/f does not depend much on the power law 
used and depends mainly on the value of the colour gradient. 
Values of r e .B used in this paper are derived from eqn. |§| 

Once the major and minor axis effective radii in both 
bands are known, we can determine the value of I e ,B from 
the K-haiid surface brightness at r e] major,B and the B — K 
colour at r etm inor,B on the dust-free minor axis. This is valid 
since B — K colour gradients on the minor axis are small 
(PB97). To convert from I e< B to the mean surface brightness 
within an effective radius we have used the relation: 



< Ie >= 10 



[0.030-(log[n]) 2 +0.441-log(n) + 1.079] 



for an r 1 /™ law (daon, Cappacioli & D'Onofrio 1993) 



2.4 Presentation of the Data 



In Table [l] we present the data used in the analysis of the FP. 
I n column (1 ) we give the New General Catalogue number 
( |Dreyer 1988; ). Column (2) lists the recession velo city of each 



galaxy in km/s, corre cted to the Local Group (Karachent- 



sev & Makarov 1996 ) , from optical heliocentric velocities 



given in the RC3. Column (3) corresponds to the scale in 
kpc/arcsec. Column (4) gives the K-b&nd geometric effective 
radius r e = ^fe^minor ■ r e , major in arcsec. Column (5) gives 
the K-h&ud central surface brightness in mag/arcsec 2 within 
r e , also from APB95, the photometric errors were assumed 
to be 0.05 (PB97). Column (6) gives the B-band effective 
radius, in arcsec, from eqn. 0. Column (7) is the B-band 
mean central surface brightness within r e in mag/arcsec 2 , 
the photometric errors were assumed to be 0.08 (PB97). 
Column (8) is the logarithm of the aperture corrected cen- 
tral velocity dispersion in km/s and its error (see Section 
2.2.2). Column (9) is the S/N per pixel of the spectrum 
used to determine the central velocity dispersion. Column 
(10) shows the parameter n of the Sersic profile and its error 
from APB95. Column (11) is the morphological type from 
the RC3. Column (12) is the disk ellipticity of the galaxy 
as given by APB95. Finally we tabulate in column (13) the 
aperture corrected ce ntral Mg2,o index in the Lick System 
(Worthey et al. 1994), in magnitudes. 



3 BULGES ON THE FUNDAMENTAL PLANE 

We determine the FP by means of an orthogonal fit to the 
function, 
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Table 2. The Fundamental Plane coefficients 



Sample Band a 



nils 



J96(Coma) B 



SrSS 0721 — -8.84 103 0.185 



Common 
fit 



B 



(0.09) (0.02) 

1.20 0.34 -9.08 28 

(0.12) (0.02) 

1.15 0.31 -8.37 46 

(0.15) (0.02) 



0.168 
0.061 

0.150 0.061 
0.098 

0.136 0.060 



P98(Coma) 


K 


1.38 


0.32 


-7.85 


64 


0.083 








(0.20) 


(0.03) 






0.128 


0.060 


Common 


K 


1.30 


0.30 


-7.31 


82 


0.092 




fit 




(0.15) 


(0.01) 






0.125 


0.063 



log(r G ) = a ■ log(cro) + P- < M +7 



(5) 



We minimise the sum of the absolute deviations of 
the points or thogonally to the rela tion, using the program 
GAUSSFIT fljefferys et al. 1987| ). We have chosen this 
method since it is relatively insensitive to the presence of 
outliers (J96, P98) and is more robust than minimizing the 
sum of the squares of the absolute deviations of the points. 

Because of the reduced size of our bulges sample, we are 
not obtaining a FP solution for our objects alone. Rather, we 
perform fits to various combinations of our sample and other 
samples in the literature, namely those of J96 in the _B-band, 
and that of P98 in the K band. In order to avoid errors due 
to distance uncertainties between clusters, we take only a 
subsample of the galaxies of J96 and P98: the Es and SOs in 
the Coma cluster. FP differences when including our objects 
provide then a measure of the deviation of our bulges with 
respect to the FP of ellipticals and lenticulars. We minimise 
any dependency on the fitting method by performing our 
own fits to the data published by B92, J96 and P98, rather 
than using their FP solutions. 

The results of the fits for B and K band data are given 
in Table |S| We list each parameter, its error, the number of 
objects employed, the rms of the solution and the rms of our 
bulges in each fit. We have also included in the last column 
the dispersion expected from observational errors; the latter 
is computed as described in Section 3.1. Uncertainties in the 
coefficients have been calculated by performing a bootstrap 
analysis with 1000 iterations. 

For the B band, the fit corresponding to the B92 sample 
excludes the dwarf spheroidals. "Common Fit" denotes the 
fit to the 28 galaxies of the Coma cluster from J96 [hereafter 
J96(Coma)] and our sample. Our fitted parameters agree 
well with the ones given by J96 , with slightly higher uncer- 
tainties due to to the smaller number of points used. 

The Coma cluster data from P98 [hereafter P98(Coma)] 
have been used to study the FP in the if-band. We per- 
formed a new FP fit to all 64 galaxies in P98, instead of the 
60 used by the authors, given the lack of information about 
which ones had been excluded by P98. Our P98(Coma) FP 
is identical to that given by P98, including the errors in the 
fit. The "Common fits" to the FP have the added value of 
extending the FP of galactic spheroids to low log(r e ) objects. 



3.1 Error Determination 

The FP fits are shown in figures |l|and|. As usual, the errors 
in log(r e ) and < fi > e are highly anticor related, with a linear 
correlation coeffi cient r of about 0.97 (|j0rgensen, Franx & 



Kjffirgaard 1995), which implies that the error in the combi- 
nation that enters the Fundamental Plane log(r e )—fl < (x > e 
are much smaller than the individual errors in log(r e ) and 
< fj, > e - For this reason those errors are not plotted. The 
error in the combination log(r e ) — /3 < fi > e was de- 
rived from a comparison with other authors (A. Graham, 
private communication). This gives us a mean value of 
cr io g (r )-^<M> =0.036. The same value has been applied for 
both bands. The observational error on the FP, including 
the uncertainties in log(cr) and distance is then 



2 



(log(r e )-/3<mu> e ) + (acriogOo)) +°d) 
and the intrinsic scatter can be determined as 



dint — [a fit 



2 il/2 
a obs\ 



(6) 



(7) 



Observational errors are listed in column 8 of Table 2. 

An independent estimate of the measurement errors 
comes from comparing our data with the literature. Unfor- 
tunately we have no objects in common with P98(Coma) 
and J96(Coma), and only one object in common with B92: 
NGC 7332. For this object, while the cr's are very simi- 
lar, we measure log(r e /kpc)=— 0.26±0.042 whereas B92 give 
log(r e /kpc)=0.010. For the surface brightness < /i > e ,s the 
values are 16.91 and 18.30 respectively. The difference in 
log(r e ) — /3 < jj, > e between these measurements, 0.70, 
is much larger than the measurement errors in either our 
or B92's data. Since the effective radius of the bulge of 
NGC 7332 is small, one might think that a difference in see- 
ing might have caused the difference. However, our sample 
contains several smaller bulges, which are lying very close to 
the FP, indicating that differences in seeing will not move 
objects away from the FP. 

We decided not to consider NGC 7332 in any of the 
common fits since it seems to be a very peculiar object and 
its observational errors are not big enough to explain such 
deviation from the different FP relations. It hasn't been in- 
cluded either in the studies of residuals of the FP with other 
parameters presented in following sections. This galaxy will 
be addressed again using, among other techniques, integral 
field spectroscopy in future papers. 



3.2 The B-band Fundamental Plane 

The distribution of galaxies in B-band FP space is shown 
in Fig. |l| Two projections are given. Fig. |l]a shows the view 
orthogonal to the J96(Coma) FP solution. Fig. jjjb shows 
the view orthogonal to the Common Fit FP (J96(Coma) 
plus this work minus NGC 7332). Substantial differences 
are found between the two distributions. Figure |l|(a) shows 
that our bulges fall slightly below the relation defined by 
the Es and SOs of J96(Coma). Other classes of objects also 
show global displacements with respect to the FP, result- 
ing in the apparent thickening of the galaxy distribution as 
we move on to smaller sizes. To analyse the behaviour of 
the different classes of galaxies we have computed the mean 
separation from the FP for each of the object groups, by 
measuring the mean deviation x along the log(r e ) direction 
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Figure 1. The B-Band FP. Velocity dispersions are in km/s, and surface brightness is in B mag/arccsec 2 . We give two different fits 
to the FP. In Fig. hp, we have performed a fit to J96(Coma) galaxies and we have also ploted our bulges and B92 galaxies following 
J96(Coma) definition. Fig. [l|b represents a common fit using all the points (J96(Coma) + our bulges - NGC7332). 



(a; = object — FP), as well as the error in the mean (Ai). 
Deviations from the FP are listed in Table ^ for each of the 
groups in Fig. [l], and for the two FP solutions (J96(Coma), 
and Common). Aperture corrections are essential to the de- 
termination of the offset. Earlier determinations by us of 
the FP offset of bulges, which did not include aperture cor- 
rections for a. yielded an insignificant offset of bulges wrt. 
J96(Coma) flBalcells 200^ ). 

Table |^ shows that our bulges show a 2-sigma devia- 
tion with respect to the J96(Coma) FP. The lower r e values 
for our bulge sample require a significant extrapolation of 
the J96(Coma) FP. The offset from the Common fit is in- 
significant, as expected since the bulge data are included in 
this fit. Interestingly, the bulges of B92 show insignificant 
deviations with respect to the J96(Coma), in apparent con- 
tradiction with B92's conclusion that bulges preferentially 
lie below the FP. We can trace the origin of this difference 
to the fact that the FP of B92, while very similar to the fits 
shown in this paper, is slightly above the J96(Coma) FP. 

The decoupling of the bright dwarf ellipticals from the 
rest of galaxies (a 15a result) is remarkable. This behaviour 
was noticed by B92 who referred to the degree of anisotropy 
as the trigger of the observed offset. The structural param- 
eters of these objects are probably intermediate between 
larger ellipticals and dwarf spheroidals. The latter are known 
to b e st ructurally different from ellipticals or bulges ( Mateo 
199S0 

We have computed the intrinsic scatter of our sample 
around the Common fit, along the log(r e ) direction using 



Table 3. Deviations from the B-band FP 



Group 


J96(C 


oma) 


Common fit 




X 


Ax 


X 


Ai 


Our Bulges 


0.074 


0.033 


-0.008 


0.032 


B92 Bulges 


0.015 


0.049 


-0.026 


0.045 


B92 Es 


-0.049 


0.018 


-0.060 


0.018 


B92 DwfEs 


-0.413 


0.027 


-0.406 


0.026 


B92 DwfSph 


-1.423 


0.217 


-1.340 


0.206 


J96(Coma) 


-0.003 


0.012 


-0.013 


0.012 



NOTES: x is the mean deviation in the log(r e ) direction (x = 
object — FP). A.x is the error in the mean. 

equation [?]. We find Oi nt ~ 0.122, which is larger than what 
can be accounted for by observational errors (0.060); see 
Table 2). This supports the idea of the existence of intrinsic 
scatter, not just for ellipticals, but also for bulges. We discuss 
clues on the origin of the intrinsic scatter in Section 5.1. 



3.3 The if-Band Fundamental Plane 

We have obtained the position of bulges on the jf-band FP 
following the same procedures used for the B-band FP. Our 
reference elliptical and SO sample here is the Coma sam- 
ple of P98 (P98(Coma)). The distribution of galaxies in FP 
space is shown in Fig. H. As for the B-band FP, two projec- 
tions are shown, one orthogonal to the P98(Coma) elliptical 
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Figure 2. The K-Band FP. 2 different fits to the FP. In Fig. |2p we have performed a fit to P98(Coma) galaxies and we have also plotted 
our bulges. Fig. bp represents a common fit using all the points (P98(Coma) + our bulges - NGC7332.). 



FP (Fig. |2p,), and one orthogonal to the Common fit FP 
(P98(Coma) plus our bulges minus NGC 7332; Fig. |b). 

Figure [2] shows that our bulges define a slightly dis- 
placed sequence with respect to the Jf-band FP of Coma 
ellipticals, extending the FP almost one decade in ra- 
dius toward low-mass spheroids. The most deviant point is 
NGC 7332. Its offset with respect to the FP is unexplained 
by measurement errors. 

The FP offsets of our bulges with respect to the 
P98(Com a) a nd the Common FP's, computed as described 
in Section 3.2, are given in Table ^. The deviations again are 
significant when using the P98(Coma) FP solution. These 
results highlight the displacement of bulges to cluster ellip- 
ticals in the K band. 

Again, the scatter of our sample around the Common 
fit in the if -Band, along the log(r e ) direction turns out to 
be larger (0.125) than the dispersion of our bulges on the FP 
due to uncer tain ties in the observational parameters (0.063) 
(see Section 3.1), implying that we find an intrinsic scatter 
in the Tf-band of 0.108. The fact that our intrinsic scatter 
in both bands is almost equal implies that for our sample 
stellar population differences are not responsible for most of 
the scatter, since this would lead to much larger scatter in 
B. We would like to remark that there is a strong paucity 
of FP data for low velocity dispersion galaxies (a < 100 
km/s), also ellipticals and SOs, and especially in the field, 
a likely consequence of the relatively recent arrival of large- 
format IR detectors. Measuring the intrinsic scatter of the 
FP in the _K"-band as well as better constraining any offsets 



Table 4. Deviations from the iC-band FP 



Group 


P98(Coma) 


Common fit 




x Ax 


x Ax 


Our Bulges 


0.054 0.029 


-0.006 0.029 


P98(Coma) 


0.002 0.010 


-0.007 0.010 



NOTES: The mean separation (x) has been determined by mea- 
suring the mean deviation x along the log(r e ) direction (x = 
object — FP). Ax is the error in the mean. 

of bulges with respect to the FP of ellipticals will require 
data for significantly larger samples. 

In Figure ^ we have plotted the residuals of the fun- 
damental plane in the B-band against the residuals in the 
if-band for the J96(Coma) and P98(Coma) fits. Residuals 
are correlated, something which we expect from the way 
these quantities were calculated. The most extreme objects 
are NGC 7332 and NGC 5879. A third outlier is NGC 5719. 
If a deviation from the plane is caused by dust, the object in 
this diagram will lie on the left of the origin. Based on pre- 
vious work, we don't think that dust extinction is affecting 
either of these 2 galaxies. The reason is the way that our fun- 
damental plane parameters have been determined. Instead 
of measuring the effective radius and the surface brightness 
in B we used the values in K and corrected them using the 
B — K colour at r e ,A' and a constant colour gradient (see 



Section 2.3.2), determined in the region where we decided, 
based on the flatness of the colour profile, that extinction 
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Figure 3. Residuals from the Fundamental Plane in B vs. those 
in K. 



was not important. In PB97 this region was determined to 
be the region outside 3" for NGC 5719 and outside 1.5" in 
NGC 5879. This method should in principle work as long as 
r e ,K would lie outside the dusty region. In fact, r e ,K = 2.77" 
for NGC 5719 and 1.90" for NGC 5879, both at the border 
of the dustfree regions, showing that the FP parameters for 
these 2 galaxies are barely affected. If the deviation is due 
to young stellar populations, the object would lie right of 
the origin. The fact that there are no objects situated here 
shows that the sample does not contain object within sig- 
nificant amounts of young stars in the bulge. The fact that 
for NGC 7332 the deviations are so similar shows that the 
deviations might be related to the details of the stellar dis- 
tribution, non-homology, or other factors. For NGC 5879, a 
small galaxy, the error in its velocity dispe rsion is so large 
that little can be concluded. (Section 2.2.2). 



4 THE Mg 2 ,o - log((7o) RELATION 
It is well known that the absorption li ne index Mg2.o corre- 



lates strongly with velocity dispe rsion (Terlevich et al. 1981 



Bender, Burstein & Faber 1993). This relation can be in- 



terpreted as a mass-metallicity dependence (e.g. Schweizer 
& Seitzer 1992, Prugniel & Simien 1996). Since, according 
to them, in early-type systems deviations are generally at- 
tributed to the presence of young stars, this relation can 
identify young stellar populations as producers of the scat- 
ter around the FP. We have determined this relation for our 
sample of bulges as well as for a number of samples in the 
literature (Fig. ^ . 

The thick solid line corresponds to a least squares fit 
to our points, while the thin solid line is a standard least 
squares fit to the bulges of B92 (solid squares). The dashed 
line is the fit given by J96 for Es and SOs. The coefficients 
of the fits are given in Table |E| We have also included the 
data of Prugniel, Maubon & Simien (2001, hereafter P01), 
a sample of 89 bulges, mostly of spiral i.e. non-SO galaxies. 
These bulges can be split into 2 subsamples: open trian- 
gles are bulges with EW(OIII)>l, filled triangles are bulges 
with EW(OIII)<l. Finally we have also added the bulges 



from Jablonka, Martin & Arimoto (1996, hereafter JMA96) 
(asteriscs). 

We first try to explain the difference between the posi- 
tion of our bulges and the ellipticals and SOs of J96. In the 
picture that the Mg2 - a relation describes the dependence 
of metal enrichment on the depth of the gravitational poten- 
tial, the position of bulges above the relation of ellipticals 
may be taken as a horizontal offset to the left - a% underes- 
timating the binding energy by not including the rotation. 
The contribution of rota tion to log (op ) may be es timated as 
0.5 log(l + 0.62V 2 /crg) ( |Prugniel fc Simien 1994"| ). An offset 
of 0.014 in Mg2 can be brought back onto the J96 relation 
by an offset of 0.071 in log(ao). A mean V/ct of 0.5 suffices 
to bring the points back on the Mg2 - a relation. We con- 
clude that non-inclusion of the rotation alone could explain 
the offset of our bulges sample from the Mg-sigma relation 
of old ellipticals. 

Fig. |i| also shows that there are significant differences 
between the various samples of bulges. To establish the na- 
ture of them we have compared the observables for objects 
in common with these authors (see Table ^j). For Mg2 the 
average differences are 0.010 ± 0.006 (RMS) for B93, -0.013 
± 0.013 for T98 and -0.024 ± 0.006 for Golev et al. The 
first 2 are consistent with the systematic error of the Lick 
system of 0.008 (Worthey et al. 1994). We conclude that the 
Mg2 values of P01 might be underestimated, leading to an 
offset in Fig. ^| One should note that, although the data of 
Golev et al. (1999) and Simien & Prugniel (1997), used in 
Table 6, were calculated in the large Lick-aperture, the data 
of Fig. 4 have been corrected to the much smaller physical 
aperture of 0.4 kpc, so that the real difference between our 
sample and P01 in Fig. 4 is probably larger. Subtle differ- 
ences in sample selection may contribute to the differences 
in Mg2 — a, despite the similar Hubble-type distributions 
and the ranges in a in our and POl's bulges. Our sample 
selection did exclude very dusty bulges (Balcells & Peletier 
1994, Peletier et al. 1999). This might have biased our sam- 
ple toward bulges without ongoing star formation. On the 
other hand, our sample is less affected by dust optical depth 
effects on Mg2 and on a, as well as from disk contamination, 
which probably affected the measurements of P01. Thus our 
sample allows more accurate measurements of stellar popu- 
lation diagnostics, but at the price of potentially excluding 
the youngest objects. We note nevertheless that our two Sbc 
bulges, at the low-<r end of the distribution, fall right on the 
relation of old ellipticals; these objects are moderately dusty. 
For log(cr) the average differences are 0.066 ± 0.013 (RMS) 
for B92, 0.084 ± 0.074 for T98 and 0.070 ± 0.025 for Golev et 
al. It appears that our velocity dispersions are consistently 
lower that the values in the literature. We think this might 
be due to the fact that our dispersions were calculated using 
minor axis spectra, while all others used major axis spectra. 
This means that in our case we did not suffer from possible 
rotation, which would tend to broaden the spectra. The shift 
of 0.07 in the log(a) direction proposed above to bring the 
bulges onto the J96 relation would bring our datapoints to- 
ward those of P01, but a much larger shift would be required 
to explain the whole difference. The data of Jablonka et al. 
(1996) fall in between us and P01, and could be considered 
consistent with us (given the uncertainties in the zero point 
of the Lick system and the above mentioned shift of 0.07 in 
the log(a) direction). 
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Table 5. The Mg 2 ,o - log(cro) Relation 



Group 


a 


A a 


b 


Ab 


J96 


0.196 


0.016 


0.257 




B92 Bulges 


0.247 


0.045 


0.258 


0.008 


Our Bulges 


0.237 


0.030 


0.271 


0.006 


P01 (OIII>l) 


0.447 


0.060 


0.183 


0.007 


P01 (OIIK1) 


0.238 


0.033 


0.201 


0.006 



NOTES: Different fitting coefficients of the Mg2,o - log(o"o) rela- 
tion, defined as Mg2,o=a ■ [log(cro) — 2.1]+b 



Table 6. Comparison with other authors 



Galaxy 


Parameter 


Us 


B93 


T98 


G99/SP97 


NGC 5422 


Mg 2 


0.307 


0.323 






NGC 5854 


Mg 2 


0.191 




0.160 


0.174 


NGC 7332 


Mg 2 


0.259 


0.263 


0.253 


0.225 


NGC 7457 


Mg 2 


0.191 




0.189 


0.171 


NGC 5422 


log(cro) 


2.205 


2.258 






NGC 5854 


log(ffo) 


1.986 




2.170 


2.017 


NGC 7332 


log(o-o) 


2.051 


2.130 


2.114 


2.117 


NGC 7457 


log(ffo) 


1.749 




1.756 


1.863 



NOTES: B93: Bender, Burstein, & Faber 1993. 

T98: Trager et al. 1998. 

G99: Golev et al. 1999 [for Mg 2 ] 

SP97: Simien & Prugniel 1997 [for log(cr )] 

Another important point to see from Fig. ^ is that the 
slope of all three bulge samples is steeper than that of J96 for 
Es and SOs. The traditional way to explain this is that low 
velocity dispersion galaxies contain a larger fraction of young 
stellar populations (P01). This can be confirmed if, e.g., the 
galaxies furthest away from the FP of ellipticals have the 
latest morphological type. For our sample the answer is not 
clear. While the 2 galaxies with type 4 are situated the most 
below the FP, no trend is visible for the others. For the 
bulges of P01, the Mg2,o deficit with respect to ellipticals is 
very pronounced, especially so for galaxies with strong [OIII] 
emission, indicating significant amounts of ionised gas. If one 
assumes that ionised gas traces ongoing star formation one 
can understand their position in the Mg2 — a relation as due 
to young stars diluting the Mg2 index. For this reason P01 
called the Mg2,o - a relation a mass-metallicity relation for 
bright ellipticals, becoming gradually a mass-age relation as 
one goes to fainter objects. Under this view of the Mg2 — a 
relation, Figure ^ is interpreted as indicating that the young- 
star contents of our bulges is much lower than in the P01 
sample, and comparable to that of SO's of B92. If indeed low 
Mg2 values are due to dilution by a young stellar continuum, 
the Mg 2 — a distribution for our bulges argues against a 
mass-age relation for small objects. 



5 DISCUSSION 

In the past bulges of spirals have been difficult galactic 
components to study. This has been due to the presence 



of dust, the accompanying disk, and the fact they are gen- 
erally small, so that very often ground-based observations 
are affected by lack of resolution. The appearance of the 
Hubble Space Telescope in th e 199 0s has vastly improved 
the situation. Peletier et al. (1999) performed a detailed 
study of galactic bulges using NICMOS observations, giv- 
ing as a result a much more detailed insight of ages and 
dust in these objects than was known before. The main con- 
clusions from this work are that early-type bulges (SO-Sb) 
behave like ellipticals and seem to be as old as cluster Es 
with very little spread in age (~ 2 Gyr). This supports the 
hypothesis of an early formation for bulges and ellipticals, 
via either dissipational collapse or early mergers. Late-type 
bulges (Sbc onwards) show hints of youth with respect to the 
early-types and could belong to a class of objects formed or 
transformed by secular disk evolution. In general they are 
smaller, younger and have lower central surface brightness. 
The question then is whether they fall on the same Funda- 
mental Plane decribed by Es and SOs. 

We have measured the mean separation of our bulges 
with respect to the FP of Es and SOs (as defined by the 
Es/SOs of J96(Coma)) and found a small but non-negligible 
offset of 0.074 ± 0.033 in the B-band and 0.054 ± 0.029 
in the Jf-band. In agreement with other authors (J96,P98) 
we find that the FP has intrinsic scatter, not explained by 
observational uncertainties. 

The strong proximity of bulges to the FP of elli ptical s 
and SOs reinforces the conclusions of Peletier et al. (199!:) 
that early-to- intermediate bulges are homogeneously old and 
structurally similar to elliptical galaxies. It is useful in any 
event to investigate the origin of the intrinsic scatter and 
the offset. We first analyse the intrinsic scatter by showing 
the FP residuals with respect to various galaxy observables. 
Then, we discuss the roles of populations, kinematics and 
homology in creating the offset with respect to Coma ellip- 
ticals. 



5.1 Analysis of the Residuals 

The existence of this intrinsic scatter has been a matter of 
debate for the last decade. For ellipticals, Mobasher et al. 
have proposed the population of AGB stars as one of 



( [1999] ) 



the sources of the intrinsic scatter of the FP in the K-band 
Structural and dynamical non-homology may be responsi- 



ble, at least in part, for the intrinsic scatter (Graham & 



Colless 1997). Recent studies suggest that the intrinsic scat- 
ter is probably a combination of age/metallicity variations 
and dynamical deviations from homology (P98 and refer- 
ences therein). For bulges differences in stellar populations 
might account for even larger scatter (e.g. Prugniel et al. 
2001). 

To obtain clues on the origin of the intrinsic scatter of 
our bulges in the FP, we analyse correlations between the 
residuals of the FP, defined as AFP = log(r e )-fit, using 
the Common fits, and different observables for our sample 
of bulges: AMj2, Galaxy Type, inclination of the galaxy 
and Sersic index n. These correlations should give us ideas 
whether the distribution of our bulges along the sequence 
described by Es and SOs depends on the presence of young 
stars [AFP-AMJ2], homology or light profile shape [AFP 
- Galaxy Type, AFP - Index of the Sersic profile (n)] or the 
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Figure 4. The Mg2,o - l°g(°"o) Relation. Filled circles correpond to our sample of bulges while solid squares represent the sample of 
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Jablonka, Martin & Arimoto 1996. Errorbars of log(cro) for our sample are also plotted. 



degree of anisotropy of the object [AFP - Inclination of the 
galaxy (i)]. 

We perform this study for the 2 bands using the Spear- 
man rank order correlation coefficients in each case. We find 
that there is no correlation between the residuals of the FP 
and any of the 3 FP parameters in either B or K, as can be 
expected from Section |j| Figure ^ shows the distributions of 
FP residuals as a function of global galaxy parameters in the 
B and jf-band. Points are labeled with the morphological 
type of the galaxy T. The probability of a null correlation is 
given within parenthesis in each panel, next to the correla- 
tion coefficient. The point labeled "-2" near the top of each 
panel corresponds to NGC 7332. This point is not included 
in the computation of correlation coefficients (see Section 

No signs of correlations are found of AFP and AMg2 
(96% and 31%, in B and K band respectively). The same 
is true if the deviations are taken from the fit of J96. If the 
deviations, both to the FP and the Mg2 - <j relation, are 
due to young stellar populations, or even to non-inclusion of 
rotational support, a correlation is expected (POT). However, 
for our sample the deviations from both relations are small, 



and this fact combined with the relatively small size of the 
sample makes us not see the correlation. 

Weak trends are found with galaxy type (64% and 12%) 
in both bands. While weak, the distribution of galaxy types 
(within our bulges) gives a suggestion that the FP residu- 
als may depend on galaxy type, in the sense that late type 
bulges tend to be below the FP, and early types above. 

We don't find a relation between AFP and n, the Sersic 
inde x. The same conclusion was reached by Graham & Col- 
less (1997) for Virgo cluster galaxies studied in the V band 
at high S/N. Such a correlation, however, could be hid den, in 
part, by stellar populations effects. Prugniel & Simien (1997) 
found a correlation between the residuals of the FP and n af- 
ter removing the scatter indiced by the stellar populations. 
Here again think the small size of the sample prevents us 
from seeing this correlation. 

Next we analyse possible correlations between AFP 
and the inclination of the galaxy (i). Projection effects and 
the presence of disks are both likely so urces of scatter in the 
FP (Saglia, Bender & Dressier 199S). Despite the limited 



range in ellipticities, we find that the probability for a sig- 
nificant correlation is 70% and 95%, in the B and K bands 
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Figure 5. Correlations of the AFP = log(r e ) — fit with global galaxy parameters. Points are labeled with the morphological type of 
the galaxy T. The number at the top left of each panel is the Spearman rank correlation coefficient. The number within parenthesis is 
the probability of a null correlation. 



respectively. J96 found that galaxies with a large ellipticity 
(presumably edge-ons) have a slight tendency to lie below 
the FP. The residual, for e = 0.6 - 0.8 » 0.10 in log(r e ), 
similar to ours. The AFP dependence on inclination can 
contribute to the global offset of our inclined s amp le w.r.t 
the J96(Coma) FP. We address this in Section 



5.2 Offsets with respect to ellipticals 



ampli 

if 



Population age/metallicity differences, internal kinematics 
and non-homology each can in principle shift bulges with 
respect to the FP of old ellipticals. We discuss each in turn. 

The position of bulges below the FP of ellipticals could 
have been produce d if bulges are on average younger than 



elliptical gal axies (Schweizer & Seitzer 1992; Prugniel & 



Simien 1996). Taking the models of Vazdekis et al. (1996) 
one finds that for solar metallicity a decrease in age from 
12.5 to 10.0 Gyr gives a shift in the log(r e ) direction of 0.06 
in B and 0.04 in K, consistent with the results presented 
here. This means that possibly the offset could be explained 
if our sample on the average would be about 2.5 Gyr younger 
than the Coma cluster samples used by J96 and P98. 

The effect of rotational sup port, although small, has 



been detected by other autho rs (Bender, Burstein fc Fabei 



199S JPrugniel fc Simien 1994 ) . Following Prugniel & Simien 



the FP via a 'rotational support term' S to be added to 
log(<7o). Using our definition of FP and after accounting for 
the inclination range for our sample, we approximate S as 



S = a 0.5 log(l + 0.62 



V' 2 

' ST, 



(8) 



with a is from Eqn. ^. Since our spectra were obtained along 
the minor axis of the galaxies, we have not been able to 
measure 5* for our sample. We can nevertheless estimate the 
amount of rotation needed to generate the observed offset. 
Non-inclusion of the rotation term offsets points downward 
on the FP diagram. A sample average of Vmax/uo = 0.66 
produces the observed FP offset. Within this model, the 
observed offset of our bulges sample with respect to Coma 
ellipticals can be explained by the non-inclusion of rotation 
data for our bulges. 

The effect of non-homology, i.e. that bulges do not all 
follow the i? 1 ' 4 law, has not been found in our data, given 
that AFP does not correlate with n, the Sersic index. As 
noted in Section 



5.1 



Prugniel & Simien (1997) did find such 
a correlation once the scatter by population effects had been 
taken into account. 



(U_994j), the rotational kinetic energy may be introduced in 



Finally, in Section 5.1 we have identified galaxy incli- 
nation as a potential origin for FP offsets. J96 model the 
inclination effect by considering the effect of rotation on the 



© 0000 RAS, MNRAS 000, 000-000 



12 J. Falcon-Barroso, R.F. Peletier, M. Balcells 



measurement of the central velocity dispersion. This effect is 
independent of the rotation support term described above. 

The previous analysis suggests that any combination 
of age and kinematics contribution to the FP can explain 
the offsets of our sample with respect to Coma ellipticals, 
whereas we fail to identify the effects of non-homology. Most 
plausibly, not a single one of the two processes is responsible 
for the entire offset. The latter probably results from a com- 
bination of the two. Nevertheless, we note that, while the 
population contribution rests on a model determination of 
ages, the kinematic offset is sure to be there; indeed, bulges 
rotate and our sample is selected to be inclined. 



6 CONCLUSIONS 

We have analysed the position of bulges of early-type spi- 
rals (type SO-Sbc) on the Fundamental Plane of early-type 
galaxies. A few conclusions can be drawn from our study: 

• We find that our sample of bulges early-type spirals lies 
slightly below the FP, as defined by E and SO galaxies, both 
in the B and K bands. The fact that bulges lie so close to 
the FP of ellipticals and SOs implies that their formation 
epoch must have been similar to that of cluster Es and SOs. 

• Both age and the unaccounted rotation effects are likely 
to contribute to the offset observed. Any single one of these 
effects can account for the observed offset, by assuming ei- 
ther an age on average 2.5 Gyr younger than the cluster ellip- 
ticals and SOs, or by assuming an average V max /ao ~ 0.66. 
Because rotation must contribute to the FP offset, the age 
differential between bulges and cluster ellipticals must be 
significantly smaller than the mentioned 2.5 Gyr. 

• We find that there is a hint that bulges of later morpho- 
logical type are situated below the other bulges in our sam- 
ple, indicating that they have slightly younger ages, again 
consistent with our results from HST colours. 

• We confirm the result of Prugniel et al. (200f) that 
bulges have a steeper Mg2,o - log(<r) relation than ellipticals 
or SOs. Contrary to current lore, the Mg2,o - a is not a 
mass-metallicity relation only, at least for bulges: Younger 
objects seem to increase the slope of the relation (POf). For 
our sample, however, the slope is shallower than for POf, 
much closer to that of the ellipticals. 
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